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A generalization of driven harmonic oscillator with time-dependent mass and 



> 

frequency, by adding total time-derivative terms to the Lagrangian, is considered. 



The generalization which gives a general quadratic Hamiltonian system does not 
change the classical equation of motion. Based on the observation by Feynman and 
Hibbs, the propagators (kernels) of the systems are calculated from the classical 
action, in terms of solutions of the classical equation of motion: two homogeneous 
and one particular solutions. The kernels are then used to find wave functions which 
satisfy the Schrodinger equation. One of the wave functions is shown to be that 
of a Gaussian pure state. In every case considered, we prove that the kernel does 
not depend on the way of choosing the classical solutions, while the wave functions 
depend on the choice. The generalization which gives a rather complicated quadratic 
Hamiltonian is simply interpreted as acting an unitary transformation to the driven 
harmonic oscillator system in the Hamiltonian formulation. 
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I. Introduction 

Time-dependent quantum mechanical systems continue of great interest. In particular, 
the system described by an explicit time-dependent quadratic Hamiltonian has attracted 
considerable attention. One of the typical examples is the harmonic oscillator with time- 
dependent mass and/or frequencies. Those studies have many applications such as in 
quantum optics , Paul trap p| , and the analyses of quantum fields in curved spacetime 
0, and they are closely related to the theory of quantum dissipation |4[]. 

In the Hamiltonian formulation of a time-dependent harmonic oscillator, Lewis and 
Riesenfeld (LR) P, |6j have shown that there exists quantum mechanically invariant op- 
erator whose exact form is determined by an auxiliary function. The invariant operator 
can then be used to find exact wave functions of quantum states. During the past several 
years, this LR method has been widely used for the study of a general quadratic hamil- 
tonian systems [0, §, ||, [T0|, PH . For the harmonic oscillator with time-dependent mass 



and frequency, the wave functions and the kernel (propagator) have been found 0, |[ |TT 



Through Heisenberg-picture approach, Ji et al. have refined the derivation, and the 
wave functions are given. The Heisenberg-picture approach have then been used fT0|j , 



with successive unitary transformations H, to find exact wave functions of the general 



quadratic Hamiltonian system The auxiliary functions in the LR method were de- 
fined through differential equations related to the equation of motion. We also note that 
the LR method can be applied for more general systems W2 . 



In this article, we will study the quadratic system in the Lagrangian formulation of 
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Feynman and Hibbs JT3[. The Lagrangian we will consider is 



L = l -M{t)x 2 - l -M{t)w 2 {t)x 2 + F(t)x + ±-{M(t)a{t)x 2 ) + ^-(b(t)x) + f(t), (1) 
2 2 at at 

where w 2 (t), F(t), a(t), b(t), f(t) are arbitrary real function of t and M(t) is also arbitrary 
real but always positive. The last three terms in the right hand side (r.h.s.) of Eq.(l) 
have no effect in classical dynamics and classical trajectory x of the coordinate x will 
satisfy the equation: 

4-(Mx) + M(t)w 2 (t)x = Fit). (2) 
at 

The most general solution of this differential equation may be composed of a particular so- 
lution and two linearly independent homogeneous solutions. The corresponding quantum 
Hamiltonian may be written as 

H = im - ami + ^ + \ Mit) ^ - wf + d(t)x + ( M) ~ m) ' (3) 



where 



c(t) = w 2 + Aa 2 -2a- 2 — a, 

w M 

d[t) = 2ab-b-F. (4) 



In the sense of differentiation of operator in Ref. ||14||, the quantum equation of motion for 
the operator x is again given as 

j t (M'x) + M{t)w 2 {t)x = F(t). 

Feynman and Hibbs have shown that the coordinate dependent part of the kernel is 
determined from the classical action (Ch.3-5 of Ref.[13|). In this article, it will be shown 



that the remained part of the kernel for the system will be completely determined from the 
Schrodinger equation and the initial condition which the kernel should satisfy. In this way 



the kernel for the system of Eq.(l) will evaluated. By the method of Ref.JT^], the wave 
functions will then be evaluated from the kernel. This Lagrangian formulation has a clear 
advantage over the Hamiltonian formulation in showing how the last three terms in the 
r.h.s. of Eq.(l) which have no effect on classical dynamics affect the wave functions. The 
classical action will be evaluated in terms of the two linearly independent homogeneous 
solutions and one particular solution, and so are the kernel and the wave functions. We 
will prove that the kernel do not depend on the way of choosing the classical solutions, 
while the wave functions are not unique and depend on the choice of the classical solutions. 



By comparing with the results on the Gaussian pure states of Ref . |16| , it is suggested that 
choosing different classical solutions might amount to acting unitary transformations to 
the annihilation operator. 

In the next section, we will consider the harmonic oscillator with time-dependent mass 
and frequency, mainly to expose our method. It will also be shown that the kernel does not 
depend on the choice of two homogeneous solutions, whilst the wave functions depends 
on the choice. In Sec. Ill, the driven harmonic oscillator will be considered. In Sec. 
IV, the system of a general quadratic Lagrangian in Eq.(l) will be considered and some 
previous errors will be corrected. The general system will be shown equivalent to the 
driven harmonic oscillator through an unitary transformation. Sec. V will be devoted to 
a summary and discussions. We add an appendix to explain how to determine the time 
dependent part of kernel from the Schrodinger equation and the initial condition. 



II. The harmonic oscillator with time-dependent mass 
and frequency 

In this section, we will apply our method to the harmonic oscillator without driving force. 
For this model the Lagrangian is written as; 

L s = )-M(t)x 2 - )-M(t)w 2 (t)x 2 . (5) 

The action (integral) from time t a to time t b is written as; 

/ Ldt (6) 

Jt a 



s 



which gives the equation of motion for the classical trajectory of the model considered 

^-{Mx) + M{t)w 2 {t)x = 0. (7) 

To find a simple expression for the classical action that is the integral along classical 
trajectory, we can rewrite the action as 

S s = -Mix \\ b -- I*" xAmx) + Mw 2 x]. (8) 
2 2Jt a [ dt K 1 J y ' 

The classical action is then simply given as 

S%(a, b) = ^M(t b )x b x b - ^M(t a )x a x a , (9) 

where x a (x b ) and x a ( x b ) denote fixed end point and ^ at t — t a (t — t b ), respectively. 

Suppose that u(t) and v(t) are two linearly independent real solutions of Eq.(7), so that 
x can be written by a linear combination of them. From the two linearly independent 



solutions, one can always find v s (t), a linear combination of the two solutions which 
satisfies v s (t a ) = 0, so that u(t) and v s (t) are again two linearly independent solutions. 
One can then easily find the quantities Q s and Q defined as 

Q s = M(t)[v s (t)u(t) -u(t)v s (t)}, Q = M(t)[v(t)u(t) — u(t)v(t)] (10) 

do not depend on time. 

The x{t) with two fixed end points x a , x b can be written as 

u{t a ) u(t a ) v s (t b ) 

Making use of this expression of x, one can rewrite the classical action in Eq.(9) as 



x a*f/+w u(t a ) u(t b )v s (t a ) x\ . ,.v 8 (tb) 



2"^" a ' 1 u {t a ) ' u(t a ) f s (t 6 ) J ' 2 v 8 (t b ) 

X a X b lM( , u Mtb) U{t b )v s {t b ) s Vsttg) , / 19 x 

+ — [M(4)( ^y - " M{ta) ^M ] - (12) 

The fact that classical dynamics are deterministic implies S c i is unique, as can be explic- 
itly proved; First, one find that S% does not depend on the scaling of u(t) or v s (t) by 
multiplying constant factors. Second, the classical action in Eq.(12) is invariant under the 
substitution of u(t) by u(t) + Cv s (t) with an arbitrary constant C. These two observations 
lead us to a conclusion that do not depend on the particular choice of u(t) or v s (t) as 
long as v s (t a ) = 0. 



Then from the formula (3-51) of Ref.[T3|], the kernel can be written as; 



K s (b, a) = exp[~(S*(a, b) + D s (t a , t b ))}, (13) 

where, as shown in appendix, D(t a ,t b ) can be completely determined from the initial 
condition and Schrodinger equation. Making use of the formulas in appendix, one can 
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find the expression of K s (b, a) in terms of u and v s as 
K s (b,a) 

V v s yi b j All 

+ M(t b 



Mjt a )v s jt a ) I u(t a ) ujt b )v s jt a ) 2 

\ 2mh v s (t b ) eXP[ 2h [M[ta)[ u(t a ) + u(t a )v s (t b ) )Xa 



^ tb) xl-2x a x b M(t a ) V ^)]. (14) 



v s (t b ) v "" y v s (t b) 

Since the kernel is uniquely determined from the classical action, our argument on the 
uniqueness of the classical action implies that K s (b, a) does not depend on the way of 
choosing u(t) or v s (t). 



To find the wave functions from the kernel along the method of Ref . JL5] , we define two 
functions p(t) and z(t) as; 



p{t) = y/v?(t)+v*{t), 
u(t) - iv s (t) 



z(t) 



Pit) 



(15) 
(16) 



After a little algebra, one can find that the kernel can be written as 
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K b (b, a) 



V^V p(t a )p(t b ) 



x exp[— ( 



L 2h K p 2 (t a ) 

oo ^n+i/ 



iM{t a )^r) + ^ 



(17) 



^ 2 n n\ n l h p(t a )'" nK \ h pit b) 
where H n is the n-th order Hermite polynomial. From now on the definition of pit) is 



modified as pit) = y m 2 (£) + t> 2 (t). From the well-known fact that 

Kix b ,t b ;x a ,t a ) = ^2i/j n ix b ,t b )^ix a ,t a ), for t b > t a , 

n 

one can find the n-th order wave function: 



n a 1 

14. 



V2~^. pit) 



uit)-ivit) ]n+lc £ [ --^ +iM{ t ) m] Hn{] 



in x . 
h'pjt) 



(19) 



which satisfies the Schrodinger equation: 

dt ~ 2M(t) &r 2 + 2 X ^' 1 Uj 

The ^ does depend on the choice of two homogeneous solutions, and any set of two 
linearly independent solutions can be used to construct the wave functions which satisfy 
the Schrodinger equation of Eq.(20). 

To have a physical interpretation of the fact that different choice of {u,v} may give 
different set of wave functions {^f, n = 0, 1,2, •••}, we consider the simplest case: 
the simple harmonic oscillator where M(t) = mo and w(t) = wq. In this case, if we 
take {u, v} as {Ccos w t, Csin w t} with arbitrary non-zero constant C, then the ip^ 
reduces to the usual stationary wave functions whose ground {n = 0) state is given as 
^ = («)i/4 e -— ^- . The choice i{ UiV ) as { cos Wot: C sin wot} with (7/ 1, however, 
gives the wave functions of probability distribution pulsating as time passes. 

In general case, by defining 7 as 71 + 272 where 

n , Mp 

7i = -rr and 72 = — z— , 
hp z hp 

we can rewrite the as 

= (21)1/4 exp(i5o(t))exp[ _l 2] (21) 
7T L 

with a real function 5q of t. Therefore, ip^ is one of the wave functions of the Gaussian 



pure states extensively studied in Ref . |T6| . There, it has been shown that any Gaussian 
pure state is the eigenstate of a certain linear combination of creation and annihilation 
operator. If we choose different classical solutions, then we could have different 7. The 



studies of Ref . [|T6f| suggest that choosing different classical solutions might amount to 
acting unitary transformations to the annihilation operator of the representation system. 



III. Driven harmonic oscillator 

In this section we will consider the system described by the Lagrangian: 

L F = ^M(t)x 2 - ^M(t)w 2 (t)x 2 + F(t)x. (22) 

Let's denote the particular solution of Eq.(2) as x p (t), so that x p (t) satisfies the equation: 

— (Mx p ) + M(t)w 2 (t)x p = F(t). 
at 

Then one can rewrite the Lagrangian as 

L F = ~[M(t)(x - x p ){± - x p )] + j t [M{t)x p {x - x p )} 

where Y(t) is defined as 

Y, p {t)= f \x p {t')F{t')dt (24) 

J to ^ 

with arbitrary constant to- The classical action S F (a, b) from time t a to tb can be written 

as 

S F (a,b;x p ) = S F (a,b) - ASi(x p (t),Y Xp (t)) \l b a 

= [-M(t)(x - x p )(x - x p ) + M{t)x p {x - x p )\ \l» a . (25) 

By adding a homogeneous solution Cu(t) + Dv(t) to the particular x p (t) with arbitrary 
constants C and D, one can have a new particular solution x' (t). By rewriting the classical 
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action as 



S%{a,b) = £[\j t (Mxt) - \j t (Mxx p ) + lj t (Mx p x) + ^x p F]dt, (26) 

one can easily find that the classical action does not depend on the choice of particular 
solution. That is, 

S F (a,b) = SZMxJ+ASfafaY,,®) \t= ~S!M^x' p )+AS l {x' p {t),Y xlp {t)) |£ . (27) 

Through the same methods of previous section, one can find the end points dependent 

part of classical action: 

5*7 , v M(t a )[x a - x p (t a )] 2 u(t a ) u(t b )v s (t a ) 
S cl {a,b;x p ) = [ — + 



2 u(t a ) u(t a )v s (t b )- 

M{t b )[x b -x p {t b )fv s {t b ) 



2 Vsitb) M t ■ t 

[x a X p (to) ) {x b Xp(t b ))- 



v s (h) 

+M(t b )x p (t b )x b - M(t a )x p (t a )x a . (28) 

The kernel can be written as [|13| 

K F (a,b) =exp[^(a,b) + D F (t a ,t b ))] = exp^S^b) + D F (t a ,t b ))}. (29) 

Since the S F does not depend on the choice of the classical solutions within the given 
restriction and D F is uniquely determined from the S F , the kernel is again unique. For 
the explicit evaluation, we require the particular solution to satisfy x p (t a ) = 0. In the 
notations of appendix, B and [3 is then given as; 

b = mm, (30 

2 v s {t b y 1 ; 

(3 = -M(t b )x p (t b )^\ - M(t b )x p (t b ), (31) 
v s \t b ) 
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and the kernel is written as 



K F (b,a) 



M(t a )v s (t a ) 
\J 2nih v s (t b ) 

x expd[^M(0(-^4 + ^ufti ) + (** - ^)) 2 M(t 6 )^ 6 



2ft L a vayv M (i a ) u(t a H(*&r v ° py ° JJ y ° J v s {t b ) 

-2x a {x b -x p {t b ))M{t a )^\ 

v s {h) 

+ 2M(t b )x p (t b )x b - 2M(t a )x p (t a )x a - M(t b )^\x 2 p (t b ) 

- r - ±„(tK(t)] 2 <ft]]. (32) 

From the expression of the kernel in Eq.(32), as in the previous section, one can find 
the n-th order wave function as 

, p , . 1 ,fi x i 1 ru(t) — iv(t),„,i Ax — xJt)) 2 , Q ,p(t),, 
^(x,t) = —== — )i—==\ w W n+ 5expp — 5TT + M t ^ 



- \ f M(z)[x p (z)^\ - i r (z)fdz\\. (33) 

In Eq.(33), {"U,f} is the set of arbitrary linear independent homogeneous solutions, and 
x p is an arbitrary particular solution. The wave functions, again, depend on the way of 
choosing classical solutions and one can explicitly find that these wave functions indeed 
satisfy the Schrodinger equation: 



IV. The general quadratic system 



In this section we will consider the general quadratic system described by the Lagrangian 
of Eq.(l). As in the previous sections, one can find the end points dependent part of 
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classical action: 

6c/(a ' 6 ' Xp) " 2 [ M^) + u(t a )v s (t b y 

M(t b )[x b -x p (t b )] 2 v s (t b ) 

2 Vsitb) Mt 

v s {h) 

+M(t b )x p (t b )x b - M(t a )x p (t a )x a 

+M(t b )a(t b )x 2 b - M{t a )a{t a )x 2 a + b{t b )x b - b{t a )x a . (35) 

The only difference of S% from is the last four terms in the r.h.s. of Eq.(35). Again, 
by requiring x p (t a ) = 0, one can evaluate the kernel in terms of classical solutions; 



K G (b,a) = 



M(tg) V S (tg) 

^ 2irih v s (t b ) 

- 2x a {x b - x p {t b ))M{t a )^\ + 2M(t b )x p (t b )x b - 2M{t a )x p {t a )x a 

v s {t b ) 

+ 2M(t b )a(t b )x 2 b - 2M(t a )a(t a )x 2 a + 2b(t b )x b - 2b(t a )x a 

- M(t b ) V ^x 2 p (t b ) 

v s {t b ) 

- f\-2f{t) + ¥&[ Xp (t)v,(t) - x p (t)v s (t)f)dt]], (36) 

Jt a V s\ l ) 

whose difference from K F (b, a) is just from the above-mentioned four terms and an integral 
of /. As in the previous sections, one can prove that this kernel does not depend on the 
way of choosing classical solutions. 

The n-th order wave function tf) n can be found from the kernel as; 
;G,„ a - 1 f n N i 1 ,u{t)-iv{t)^ h 



COM) 



Vz^J. J^t) v pit) 



x exp[UM(t)a(t)x 2 + (M(t)x p (t) + b(t))x]] 
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xexp[ 2h { p*(t) +lM{t) p(ty lHn l h P (t) j 

XeXPi ^ [ 2 v(tf p{t) 

+ /'(/(*) - ^Mz) V -^\ - x p {z)f)dz\\. (37) 

Jt I V[Z) 

Again, {u, v} is the set of arbitrary linear independent homogeneous solutions, and x p 
is an arbitrary particular solution. One can explicitly apply the Schrodinger equation to 
this wave functions, to find that they indeed satisfy the equation. 

In the Lagrangian of Eq.(l), the conjugate momentum p of the coordinate x is written 
as p = Mi + 2Max + b. One may interpret x p as the classical coordinate, and the 
conjugate momentum is then written as 

p p = Mx p + 2Max p + b. (38) 

As in the Sec. II, we define a 7' as j[ + i^' 2 where 

7 =71 and 72 = — J" 2a + ^ . 

n p 

Then, the wave function tp^ can be simply written as 

t) = _^=(^)i e «(«) exp[-^(x - x p f + Up^H^ix - xp)), (39) 



where 8(t) is a real function of t. The wave functions agree with those in Ref.[|TTJ except 
the fact 5(t) is real which is necessary for the conservation of total probability J ip^*^dx. 
The expression of ip^ i n Eq.(39) shows that i/jq is a wave function of a Gaussian pure 



state |IB| , so the discussions of Sec. II are still valid in the general case. 

With the wave functions, one can calculate the expectation values of operators, and 
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the uncertainty relations read; 



,< (Ax) 2 > m m < (Ap) 2 > m = ( m < X 2 > m - m < X > 2 m )( m < p 2 > m - m<P >m) 



(m + l -fn 2 [l + ^(2Map 2 + Mpp)\ (40) 



m+ i< (Ax) 2 > 

m m +1< ( A Pf >m 



-7=(- — p^- * (m + w) 3 , ^ 1 }[2Ma + M'- + i— 

V2 K n J K \j(m+l)h(u + iv) p p 2 



x [Pp - l^-^r 1 ^ + w){2Ma +M - P + V )L (41) 



1. 

+2 < (Ax) 2 > +2 < (Ap) 2 > 

= (m + 2)(m + 1)(^) 2 (« + *t>) 4 (2Ma + M- + i^) 2 , (42) 



A)2( u + it ,)<( 2 Ma + M£ + i^ 2 
2s I P P 



with the notation n < O > m = J^ 00 ip^*(x,t)Oip^ l (x,t). If we take u = pcos9 and v = 
psinO, then il = Mp 2 # and the functions p(t), 9(t) should satisfy: 

+ 2^0+^0 = 0, 
p M 

P + M P ~ P + P = 
With these notations, the uncertainty relations in Eq. (40-42) are written as; 

m <(Ax) 2 > m m <(Ap) 2 > m =(m + i) 2 n 2 [l + -J-(2a+^) 2 ], (43) 

2 6 2 p 



m+1 < (Ax) 2 > m m+1 < (Ap) 2 > m 



4 0> l J( m +l)h 1 U P J 



x [ 2a + £ + i0 - - 2 ^ 26 Pp e- ie ] (44) 
P yJ(m + l)Mh 

. l+2 < (Ax) 2 > m m+2 < (Ap) 2 > m = {m + l) \ m + 2) % 2 e^l[2a + t + i 9] 2 , (45) 

4 9 2 p 
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respectively. The uncertainty relations of Eq. (43,45) exactly agree with those of Ref.fTT 
but the uncertainty relation of Eq.(44) differs from the corresponding one there. 

The terms which do not affect the classical dynamics of the model in Eq.(l) can be 
written as L — L F . The effects of those terms on the wave functions could simply be 
represented by writing ijj^ as; 

^{x, t) = exp[^ J\l - L F )(x, ^, z)dz] x if; F {x, t). (46) 

This relation suggests that ip^ can be obtained from ip F by acting unitary operator U ; 

U = exp[^{M{t)a{t)x 2 + b{t)x + J* f{z)dz)\. (47) 

By defining operator Op, Oq as 

F = -ih- + H F , G = -ih- + H, (48) 
at ot 

one may find the relation 

UO F U ] = G , (49) 

which proves that the Schrodinger equation of general quadratic system is equivalent to 
that of the driven harmonic oscillator through the unitary transformation. 



V. Summary 

The Feynman and Hibbs formulation (or an observation) on the quadratic Lagragian sys- 
tem gives a good explanation on the fact that the quantum wave function can be written 
in terms of solutions of classical equation of motion. By developing the observation, we 
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find the kernel and wave functions of the general quadratic system in terms of classical 
solutions. Furthermore, the kernel is shown to be independent from the choice of clas- 
sical solutions. These results are then used to show that the general quadratic system 
is equivalent to the driven harmonic oscillator through a unitary transformation. This 
fact [jnj] shows that unitary transformation (or, canonical transformation, its classical 



correspondent) could make the problem simpler or more complicated, and could change 
the uncertainty relations as in Eq. (40-42). 
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Appendix 

In this appendix, it will be shown that, if 

K(a, b) = exp[-(Axl + Bx\ + hx a x b + ax a + (3x b + s)] (A.l) 

where A, B, h, a and (3 are already known functions of t a , t b , then the function s(t a , t b ) is 
uniquely determined from the initial condition and Schrodinger equation. For the system 
of Hamiltonian given in Eq.(3), the kernel should satisfy the Schrodinger equation: 



d 1 ft d 2 h d M(t b )c(t b ) 2 

''aT K = om/i V V~^~) _ 2a{t b )-x b — + r £ ft 

ot 6 2M(t b ) i ox b i ox b 2 



6(4) U + d{h)xb + ( 6*fo) _ + - Mt6))] ^ (A _ 2) 



M(t b ) z dx b v y V 2M(4) 

which gives the following differential equations; 

OA h 2 
dt b ~ 2M{t b ) 

8B 2B 2 A , M(t b )c(t b ) 



(A.3) 



+ 4a(t 6 )£- v ^^" 7 (A.4) 
+ 2a(t 6 )fc (A.5) 
/i (A.6) 



54 Af(t 6 ) 

94 ~ ~M(t b ) 
da hp b(t b ) 

df b = ~M(t b ) + M(t b ) 

I - -ii^r^ + M/-mt) +m - ma ^ (A - 8) 

With the explicit expressions of A,B,h,a and (3, one may check that Eqs.(A.3-7) are 
satisfied. For example, in the general quadratic system considered in section IV, B is 
given as B(t a ,t b ) = + M(t b )a(t b ) which satisfies Eq.(A.4). 
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s{t a ,h) can be determined from Eq.(A.8) up to a function g(t a ,t ); 

h fa B 



IL f 

S(t a ,t b ) = g{t a ,t ) - ~ 

I Jt, 



-dt 



t0 M(t) 

~ [tTTTTT - + TTTTTT ~ /(*) + tha{t)}dt. (A.9) 

Jt a l 2M(t) M(t) 2M(t) 

A wave function ip satisfies the integral equation 

/oo 
K(x b , t b ; x a , t a )ip(x a , t a )dx a . (A.10) 

In the limit of t b — > t a , the classical actions approaches to 



2(t b - t a ) 

In order that the relation of Eq.(A.lO) be satisfied in the limit, the kernel should satisfy 
the relation: 



K(b, a) 



M(t a ) _ r iM(t aJ 



\2irih(t b -t a ) iy 2h(t b -t a ) 
This initial condition determines the g uniquely. 



ex P[7^7^ —r(xa-Xb)] as t b — > t a . (A.ll) 
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